Is the future universe singular: Dark Matter versus modified gravity? 
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The fundamental problem of the occurrence/removal of finite-time future singularity in the uni- 
verse evolution for coupled dark energy (DE) is addressed. It is demonstrated the existence of the 
(instable or local minimum) de Sitter space solution which may cure the Type II or Type IV future 
singularity for DE coupled with DM as the result of tuning the initial conditions. In case of phantom 
DE, the corresponding coupling may help to resolve the coincidence problem but not the Big Rip 
■ (Type I) singularity issue. We show that modified gravity of special form or inhomogeneous DE 

fluid may offer the universal scenario to cure the Type I, II, III or IV future singularity of coupled 
f^*) , (fluid or scalar) DE evolution. 
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I. INTRODUCTION 

The singularity issue has the fundamental importance in the modern cosmology. The early universe may appear- 
totally different what depends on the presence/absence of the initial singularity. There are various points of view on 
this issue. Generally speaking, the singularity problem may be properly understood/resolved only in full quantum 
gravity theory which does not exist so far. 

With the discovery of the late-time acceleration the singularity problem became even more important. The reason 
is that observational data favour the ACDM cosmology with the equation of state (EoS) parameter being very close to 
— 1. This means that phantom/quintessence Dark Energy (DE) models with effective EoS parameter w approximately 
equal to —1 are not excluded. Unfortunately, it was discovered [l[ that many of such DEs may lead to one of four 
different finite-time, future singularities in the universe evolution. Definitely, the presence of finite-time, future 
singularity may cause various problems/instabilities in the current black holes and stellar astrophysics. Hence, it is 
very interesting to understand if any natural scenario to cure such singularities exist. The purpose of present work 
is to consider the classical scenario of the future singularity removal for coupled DE. In the next section we consider 
, phantom DE coupled with Dark Matter(DM). It is demonstrated the existence of instable de Sitter solution which may 
solve the coincidence problem but does not cure the future singularity. Third section is devoted to the consideration 
of DE fluid which may develop all four types of future singularity. It is shown that its coupling with DM may cure 
Type II and Type IV singularities but not Big Rip (Type I) and Type III ones. In fourth section it is demonstrated 
that special models of modified gravity or of the inhomogeneous EoS fluids [2] may offer the universal scenario to cure 
the future singularity of any type. This is demonstrated for fluid or scalar-tensor singular DE. Some outlook is given 
in the conclusion. 
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II. DECAYING PHANTOM DARK ENERGY: THE SOLUTION TO THE COINCIDENCE PROBLEM 

Let us remind several simple facts about coupled phantom DE. The DE conservation law is given by 

p DE + 3H {I + w) p BE = 0, (1) 

while the first FRW equation 

~^H 2 — pde , (2) 

admits a solution ( w is less than —1) 



H = ^± , (3) 
t, - t ' W 
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which has finite-time future singularity, called Big Rip, at t = t a (for number of earlier works studying singular 
phantom DE era, see 0,0]). 

We now consider the model where phantom DE couples with dark matter. The conservation law is 

Pde + 3H (1 + w) pde = -Qpde , Pdm + 3i?p D M = Qpde • (4) 

Here Q is assumed to be a constant. The first equation can be solved as 

PDE = PDE(o)a- 3(1+w) e- Qt . (5) 

Here Pde(o) is a constant of the integration. Then the second equation Q can be solved as 

Pdm = Qa(t)- 3 J dt'p nm a- 3w e- Qt . (6) 
The second FRW equation is now given by 

' 2H + 3H 2 ) = p = wpuE = wpvE^a-^+^e-^ . (7) 



K \ 



An exact solution of is de Sitter space 



„(/i = „„<-™t* ( g = - _^_ ) (S) 



where ao is given by 



Q \ -3(l+io) , n , 

' ■= wpuE(Q)a ■ (9) 
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Note that we are considering phantom DE with w < — 1 and therefore H {J5| is positive and Eq. (|9j) has a real solution. 

Hence, the coupling of the dark matter with the phantom DE gives de Sitter solution instead of Big Rip approaching 
solution ([3]). This does not always mean that Big Rip singularity could be avoided but gives a possibility that the 
universe could evolve to de Sitter universe instead of Big Rip. 

One may identify the Hubble rate H with the present value of the Hubble rate H = — 37^^ = Hq ~ 10~ 33 eV. 
Eq. (HJ shows that DE density is a constant 

-3(1 + 10) /,A\ 

PDE = PDE(0) a ■ ( 10 ) 

Then Eq.© can be integrated as 

pDM = PDMoa - (1 + w) /0 DE (o)a (11) 

Here puu is a constant of the integration but the first FRW equation 

3 

— H 2 = p D M + PDE , (12) 

shows pdmo = and therefore the dark matter density /Odm is also constant: 

Pdm = -(l + w) Pde(o)Oo 3{1+w) =-(1+w) p DE ■ (13) 

Then if the de Sitter solution JH| is an attractor, by choosing 

1 4 
+ i.e. w~--, (14) 

the coincidence problem could be solved. That is, even if we start with a wide range of the initial conditions, the 
solution approaches to the de Sitter solution, where the ratio of DE and DM is approximately 1/3 almost independent 
from the initial condition. 

If DM does not couple with DM, the densities behave as pdm ~ a_3 and pde ~ a~ 3 ( 1+w ' , that is, DM density 
decreases but the DE density increases when the universe expands. This requires the fine-tuning for the initial 
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condition of the densities so that the density of DM has almost same order with that of DE. This is so-called 
coincidence problem. 

In order to investigate if the de Sitter solution |(8|) is an attractor or not, we consider the perturbation as 

a = a e-^ t+A «. (15) 
Here A(t) is assumed to be small. The second FRW equation J7]) gives 



1 (ok _JQ_a\ a -3(1+-) a 3(1+^0 / Q \ 



2 

A, (16) 



which is very simple linear differential equation with constant coefficient. Here Eq.© is used. By assuming A = e At , 
we find 

= A 2 - T ^-A + 3(1 + W ) ( on Q X , (17) 

that is 



1 + w \3(l + w) 




Then A_ < but A + > and therefore the de Sitter solution ([8]) is not stable. Since A_ < 0, however, the solution 
is saddle point and therefore if we choose the appropriate initial condition (not by fine-tuning), there is a solution 
approaching to the saddle point de Sitter solution ((8]) . Although we need to study the global structure of the space of 
the solutions, such the initial condition could correspond to the direction towards to A_. Hence, generally speaking, 
the coupling of DE with DM does not prohibit the existence of the (Big Rip) singular solution. Then if the singularity 
corresponds to the stable solution, since the de Sitter space is not completely stable, although there is a solution 
approaching to the de Sitter solution, the solution will finally evolve to the singular solution. Hence, the appropriate 
choice of initial conditions may help to realize the non-singular de Sitter cosmology which also solves the coincidence 
problem. At least, one may expect that the future singularity occurs at more late times than without coupled DM. 



III. COUPLING OF DARK ENERGY WITH DARK MATTER: SINGULARITY AVOIDANCE? 



Let us consider the simple example of perfect fluid with the following equation of state (EoS) @ : 

P = - P + Ap a , 



with constant A and a. We work in the spatially flat FRW space-time : ds 2 = —dt 2 
the Hubble rate is found to be 



H = { 



t -t 



V3KAt 



Be' 

Ct l/(l-2 a ) 

or C (t -t) 



l/(l-2a) 



when a = 1 , A > 
when a = 1 , A < 
when a = \ , A < 

when a^l, \ 



(19) 

(*) 2 £i=i, a ,3M a - Then 

(20) 



Here B, C, and C are positive constants. Now one can describe the future, finite-time singularities of the universe 
filled with above dark fluid for different choices of theory parameters (see ref.@). For more complicated DEs leading 
to all four types of future singularity, see @, 0|. When a < 0, there occurs Type II or sudden future singularity @, Q- 

When < a < 1/2 and 1/(1 — 2a) is not an integer, there occurs Type IV singularity. When a = 0, there is no any 
singularity. When 1/2 < a < 1 or a = 1 and A < 0, there appears Type I or Big Rip type singularity. When a > 1, 
there occurs Type III singularity (for classification of all four types of future singularity, see [l(). 

In case of Type II singularity, where a < 0, H vanishes as H ~ (t a — t) 1 ^ 1 ^ 2 ^ when t — * to and therefore we find 



p vanishes as it follows from the FRW equation 
to 



(3//v 2 ) H 2 = p. Then, near the singularity, the EoS lfT9|) is reduced 



P ~ Ap a . 



(21) 
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On the other hand, in case of Type I singularity, where 1/2 < a < 1 or a — 1 and A < 0, H and therefore p diverges 
when t — > to. Then the EoS (fT9|) reduces to 



p povp (1-A)p. (22) 

In case of Type III singularity, where a > 1, H and p diverge when t — > t and therefore the EoS lfl9|) reduces to 

p ~ Ap" . (23) 
We now consider DE (fl9l) coupled with DM as in (@]). Then the conservation law is given by 

Pde + 3#Apg E = -Qpbe , Pbm + 3HpuM = Qpbe ■ (24) 

The solution of Gl is 



PdeW = e- Qt (-3A (1 - a) J dt' H(t')e^ Qt '^ ' " , p DM = Qa(t)- 3 j dt'a{t'f p BE (t') . (25) 
Then if A < (and Q > 0) , there is a de Sitter solution H = Hq with a constant -£/q , which is given by solving 



^ = ( 1+ Jo)Hr> 



and pde and pdm are constants 



PDE= vT^rJ ' pdm = sio r-Q-J • (27) 

This demonstrates that if H ~ Q, we find Pdm/pde ~ 1/3 and the coincidence problem may be solved. If H = Q, 
Eq. ljM)) determines the value of A: 

^=-Ki>°P- < 2s » 

One now investigates the (in) stability of the de Sitter solution H = H by putting H = H + SH. The perturbation 
of energy-density is 

S PBE = -3A (-^^) 1 a e^ 1 "^' J* dt/SHffleW* . (29) 
The second FRW equation becomes 

_ _L (sH + 6ff «ff) e (1 -« )c?t ' = 3A (l + ^ ("^^) ^ / dt' SHit')^ 1 -^' . 
By differentiating the both sides of ([30| . we find 

= <5# + {6i/ + (l-a)(3}<5ij+{6 J ffo(l-a)g + 3A K 2 fl + ^-) f-^T^V " 1 SH 



(30) 



3Jj y V Q 

= SH + {6H + (1 - a) Q} SH + 3H Q (1 - 2a) QSH . (31) 
In the second equality, Eq. (f26|l is used. Assuming SH oc e At , one gets 

= A 2 + {6H Q + (1 - a) Q} A + 3iJ (1 - 2a) Q , (32) 

whose solution is given by 

A = A± = _ 6g o + (l~a)Q ± 1 | {6jfff) + (1 _ a) Q} 2 _ (1 _ 2a) Q y (33) 
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Then since Hq, Q > 0, if 



a < \ , (34) 



both of A± are real and negative if 

D = (6-ffo + (1 - a) Qf - 12H (1 - 2a) Q > , (35) 
or complex but the real part is negative if 

D = (6ff + (1 - a) Qf - 12H (1 - 2a) Q < . (36) 

Then as long as a < | ([34|) . the de Sitter solution is stable and therefore the singularity could be avoided. On the 
other hand, if 

6H + (l-a)Q >0, 3ff (l-2a)Q < 0, (37) 

that is 

6ff 1 

l + y>a>-, ( 38 ) 

we find A + > and A_ < 0, as in the case of the previous section. In fact, the previous section corresponds to a = 1 
case. If 

6iTo + (1- a) Q <0, 3£T (l-2a)Q < 0, (39) 

that is 

6-ffo , x 

a > 1 + , (40) 

we find A± > and the de Sitter solution is totally unstable. 

The existence of the singular solution itself, even in the presence of coupled dark matter, can be confirmed by 
substituting the singular solution without DM into the equations. 

Hence, Type II singularity, where a < 0, and Type IV singularity, where < a < 1/2 and 1/(1 — 2a) is not an 
integer, could be cured by the coupling of DE with DM. As we mentioned, even if there is a coupling of DM with DE, 
there could be a singular solution. The de Sitter solution with a < \ is, however, at least a local minimum. Then 
if universe started with an appropriate initial condition, the universe evolves into the de Sitter one (asymptotically 
de Sitter universe). Type I (Big Rip) singularity, where 1/2 < a < 1, and Type III singularity, where a > 1, could 
not be removed by the coupling of DM with DE. Even if a solution goes near the de Sitter point, the solution could 
evolve into the singular solution if the singular solution is stable. Nevertheless, the fact of avoidance of some future 
singularities due to coupling of DE with DM looks quite promising. 



IV. MODIFIED GRAVITY CURING THE SINGULARITY 



Since the coupling of DE with DM does not always remove the singularity, we now consider what kind of the fluid 
could cure the future singularity. In case of the Big Rip singularity, for example, the energy-density of DE diverges 
like pde ~ 1/ (t s — t) 2 ~ R when t — > t s . Here R is the scalar curvature. Then one is interesting in a fluid whose 
pressure is positive (and energy-density is positive) and grows up more rapidly than DE pressure. There is no such a 
fluid with constant EoS parameter. However, one can consider the pressure which is proportional to a power of the 
curvature, for example, 

Muid ex R 1+e , (41) 

with e > 0. Then the total EoS parameter becomes greater than —1 for large curvature and Big Rip does not occur. 

This kind of inhomogeneous effective fluid 0] could be realized by quantum effects (for instance, taking account of 
conformal anomaly) or by modified gravity (for general review, see [9j). As an example, we consider F(R) gravity, 
where F(R) = R + f(R) behaves as f(R) cx R m . 



6 



When R is large, if m > 1, the contribution from matter, DM and DE could be neglected. There occurs the 
following solution: 



H 



(m-l)(2m-l) 
m-2 



(42) 



which gives the following effective EoS parameter: 



2H 
3~fP 



2(m- 2) 



3(m- l)(2m- 1) 



(43) 



which is greater than — 1 if 2 > m > 1 or m < 1/2. In case of m < 1/2, however, the Einstein-Hilbert term i? 
and/or the dark matter could dominate over f(R) term and we do not consider this case. In case of 2 > m > 1, 
since w c g > — 1, there does not occur the Big Rip singularity or any kind of the future singularity. We should note 
that there occurs the Big Rip singularity when m > 2. In case of m = 2, as Eq. (|42|) tells, the power law solution 
like Big Rip singularity is prohibited. We can find that there appear (asymptotic) de Sitter solution, instead of the 
power law solution as in (|42|) . which is consistent with lj43"l) . and therefore the singularity is prohibited. Therefore, 
the introduction of special form of /(i?)-term prevents the future singularity. Note that as other DEs, the modified 
gravity itself may lead to all four possible future singularities [7j which may be cured b y R -term 0, Efll (for related 
discussion of curing the Type II singularity in special modified gravity by R 2 term, see 
On the other hand, if we choose 



iy it -1 

m 



m = m± = 



7±V73 



(44) 



we find w c ff vanishes, which corresponds to the dust-like dark matter or usual matter. 
Let us consider the following example: 



F(R) ~ f+R m ° + f-R r ' 



(45) 



If we choose mo to be 2 > mo > 1, the first term will prevent the singularity when the curvature is large and the 
second term might behave as dark matter when the curvature is small since m_ < 0. In (@3|), if m > 2, there could 
occur Big Rip type singularity. In order that H > in lj42j) . one may assume t < at present universe, or equivalently 
shift t as t —>■ t — to and rewrite 11421) as 



H 



(m-l)(2m-l) 



t -t 



Then if t < to in the present universe, H is positive and there occurs the Big Rip singularity at t 
hand we may consider the scalar tensor-theory whose action is given by 



(46) 

to. On the other 



dx 4 



V e 



(47) 



with constants Vq and 4>o- The action (|47| admits the following solution 



= 0o In 



H 



2^2 



4t 



1 - 



2V n 



(48) 



Then if 



2 j.2 



(49) 



we find 



scalar-tensor theory (|47j) coupled with f(R) oc R 



and the case of l|49j) corresponds to the dark matter. We may consider the model, the 
gravity with m > 2 like the Brans-Dicke theory as 



5st = J dx^y/^g 



R I 

2^ + foRm ~ 2 S " 



Vbe *o — Uoe't'a R 



(50) 
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with a coupling Uq. The term with Uq could express the interaction between gravity and the scalar filed. The action 
(|47l) admits the de Sitter solution as we see now. Assume the scalar field cf) and the curvatures is covariantly constant: 

4> = c, R = i? , R^v = jRog^u , (51) 

with constants c and Rq. Then the equations corresponding to the Einstein one and the scalar field have the following 
form: 

= ^ + (2-m)f R™-2Voe-$S - U e^ R , 
= F e"W - U e^Ro . 
Combining Eqs. (|52"ll and l(53|) . one obtains 

(54) 

/ Vn \ 2 (4m-2)c 

m)/ f e — = 3U . (55) 

If Vo, Uo > 0, /o < 0, and m > 2, the l.h.s. of ([55| is positive and monotonically decreasing function of c and the 
l.h.s. vanishes in the limit of 4> — > +°o but positively diverges in the limit of 4> — > — oo. Therefore (|55|) has a unique 
solution with respect to c. Then Eq . (|54)1 shows R is constant and positive, which expresses de Sitter universe. We 
should note that when Uq — ► 0, R diverges and therefore there is no de Sitter solution without the coupling Uq. 

Hence, we demonstrated that DE with inhomogeneous EoS corresponding to special form of modified gravity may 
easily cure the future singularity of any type. It indicates that if our universe does not like the future singularity, 
then modified gravity should play the role of DE. 

V. DISCUSSION 

In summary, we considered DE model which may lead to all four types of future singularity in the late-time universe 
evolution. It is demonstrated that its coupling with DM may cure Type II and Type IV singularities (but not Big 
Rip and Type III) already on the classical level. It turns out that only modified gravity /inhomogeneous EoS DE may 
suggest the universal classical recipe to remove any of the known future singularities. This is shown for fluid DE as 
well as for scalar DE coupled with modified gravity. Moreover, there are viable non-singular modified gravities which 
describe the unification of the early-time inflation with late-time acceleration. Adding such non-singular modified 
gravity to singular DE model effectively removes the future singularity of any type as is described in ref.[I3|. In this 
respect, the alternative gravity DE may seem to be more fundamental theory than more traditional (scalar, fluid, 
etc.) DE. 

As final remark, one should stress that our consideration is totally classical. Nevertheless, it is expected that 
quantum gravity effects may play the significant role near to singularity. It is clear that such effects may contribute 
to the singularity occurrence/removal too. Unfortunately, due to the absence of complete quantum gravity only 
preliminary estimations may be done. However, already the account of one-loop quantum gravity effects indicates to 
the possibility of the removal of future singularity This gives one more argument in favor of modified gravity as 
the possible universal regulator of future singularity. 
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